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This paper examines the classical dynamics of false vacuum regions embedded in surrounding 
regions of true vacuum, in the thin-wall limit. The dynamics of all generally relativistically allowed 
solutions - most but not all of which have been previously studied - are derived, enumerated, 
and interpreted. We comment on the relation of these solutions to possible mechanisms whereby 
inflating regions may be spawned from non-inflating ones. We then calculate the dynamics of 
first order deviations from spherical symmetry, finding that many solutions are unstable to such 
aspherical perturbations. The parameter space in which the perturbations on bound solutions 
inevitably become nonlinear is mapped. This instability has consequences for the Farhi-Guth-Guven 
mechanism for baby universe production via quantum tunneling. 

PACS numbers: 98.80.Hw, 98.80.Cq 



I. INTRODUCTION 

Nearly two decades ago, a series of papers Q, 0, Q 
began to investigate the possibility of creating an in- 
flationary universe "in a laboratory" - that is, inside a 
surrounding region of much lower vacuum energy. The 
spacetime of such a "bubble universe" was modeled as a 
spherically symmetric de Sitter region (the false vacuum) 
joined to a surrounding Schwarzschild geometry (the true 
vacuum) by an infinitcsimally thin domain wall. 

These studies found that the inflating (false vacuum) 
region could not avoid collapse unless either the null 
energy condition or cosmic censorship were violated in 
the full spacetime 0, but that the creation of an en- 
during inflating region might be possible via quantum 
tunneling 0, 0- I n this picture, henceforth denoted 
the Farhi-Guth-Guven (FGG) mechanism, a classically 
constructable expanding bubble, which would classically 
re-collapse (both from the inside and outside perspec- 
tive 0), instead tunnels to a new solution in which the 
inflating interior expands forever, while an outside ob- 
server sees a black hole (3^. The probability of this oc- 
curring can be calculated using the techniques of semi- 
classical quantum gravity, and is extraordinarily small. 

Though nearly-miraculous, this process has garnered 
new interest recently, primarily because of evidence that 
our universe may have a fundamental positive cosmolog- 
ical constant. If so, it will asymptotically approach ev- 
erlasting equilibrium as de Sitter spacetime. Given eter- 
nity, even the most unlikely process - such as the creation 
of inflating bubble universes - will eventually occur. Tak- 
ing this one step further, our observable universe could in 
fact be such a bubble universe, which arose from equilib- 
rium de Sitter space and is currently returning to it. This 
would realize the old idea of Boltzmann that the universe 
is fundamentally in equilibrium, but that extremely rare 
downward fluctuations in entropy periodically occur and 
allow the transitory existence of non-equilibrium regions 
that see entropy steadily increasing. 

The classic problem with this idea was pointed out in 



its new context by Dyson, Susskind, & Kleban [5j: if our 
observable universe resulted from a downward entropy 
fluctuation that evolved with increasing entropy to the 
present time, then it is vastly more likely for this to have 
occurred by a fluctuation to our observable universe ten 
minutes ago (replete with incoming photons and mem- 
ories in our brains to convince us that it is older) than 
by a fluctuation all the way to the much lower entropy 
corresponding to inflation |36j. Albrecht and Sorbo @], 
however, argue that inflation might avoid this problem by 
turning a tiny region of low entropy density into a very 
large one. Their calculation shows that the FGG mech- 
anism requires a much smaller entropy fluctuation than 
directly creating the large post-reheating region that will 
result from it, thus resolving the paradox. 

There may, however, be reasons to doubt that the cre- 
ation of small regions of false vacuum and subsequent 
tunneling are plausible events. First, there is no regular 
instanton describing either the nucleation of small regions 
of false vacuum or the gravitational tunneling event; 
there is such an instanton describing tunneling to false- 
vacuum, but over a huge regi on, larger than the true- 
vacuum horizon (see, e.g., [HlEI). Second, Banks 
has pointed out that it is hard to understand the tunnel- 
ing process holographically - it appears that the inflating 
region inside the bubble has many more states than the 
black hole it is "contained in" (see also 0]). Third, 
Dutta & Vachaspati Q have recently (and previously 
in argued that general causality considerations pre- 
clude the formation of a small false- vacuum region inside 
a large true- vacuum one. 

Both in terms of the initial conditions for inflation, and 
also the more general issue of what processes can lead to 
transitions between vacuum states (as is important in 
understanding the string theory "landscape"), it is cru- 
cial to understand bubble universes, whether they can 
form, and with what probabilities. To this question the 
present paper makes the following contributions: first, 
we organize and interpret all of the thin-walled, spheri- 
cally symmetric one bubble spacetimes. We then show 
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that expanding false vacuum bubbles are unstable to non- 
spherical perturbations; if these bubbles start at a suffi- 
ciently small initial radius, then they inevitably become 
nonlincarly aspherical before tunneling might occur in 
the FGG mechanism. It is unclear in this case whether 
tunneling to an inflationary universe inside the bubble 
can occur at all, or with what probability. 

The plan of the paper is as follows. In Sec. [HI the 
allowed solutions to the junction conditions are enumer- 
ated, putting the previous work into context. We pro- 
vide a concise reference for all of the possible spacetimes 
with one false- vacuum bubble and arbitrary positive cos- 
mological constant, then discuss the existing, and some 
new, interpretations of these solutions. In Sec. GUI we de- 
rive the first order perturbation equations, and demon- 
strate the existence of an instability in bound solutions. 
In Sec. |W]we integrate the equations to investigate the 
parameter space for which non-linear perturbations are 
unavoidable, and we conclude in Sec. 



II. JUNCTION CONDITIONS 

The dynamics of inflati ng regi ons has been discussed 
by a number of authors |4 lltlll7lll8| |37| using the junc- 
tion condition formalism. These works study a spheri- 
cally symmetric region of false vacuum (high energy den- 
sity) in a surrounding region of true vacuum (lower en- 
ergy density). The wall separating the regions is assumed 
to be very thin compared to the radius of the region of 
false vacuum. One can obtain the dynamics of the wall 
by requiring metric continuity across the wall and then 
solving Einstein's equations. Under spherical symmetry, 
the dynamics of the problem then reduce to those of the 
bubble wall's radius. This radius is a gauge invariant 
quantity because it simply quantifies the curvature of the 
bubble wall worldsheet, and any observer can measure 
it by comparing the normal to the wall at two nearby 
points. 



A. Interior and Exterior Spacetimes 

Let A + be the cosmological constant in the true vac- 
uum region. Then if A + > 0, the region is Schwarzschild- 
de Sitter spacetime with metric 
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in the static foliation. Fixing A + , there are then three 
qualitatively different casual structures characterized by 
the value of M (see H^), due to the nature of the three 
roots of a s ds(0- 
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FIG. 1: Conformal diagram of the Schwarzschild de Sitter 
geometry for 3M < A+. 



form: 
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r„ = 2(A+)- 1 / 2 cos^- + — 
where 

cos 6» = -3Af(A+) 1/2 , 
and 7T < 9 < 3ir/2. We can label them as 
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And the range of 6 means that they lie in the ranges 
r nGg < < 2M < r BH < 3M < r c . 

The two positive roots correspond to the black hole 
and cosmological horizons. The conformal diagram for 
this spacetime is shown in Fig. (See [2(j for a demon- 
stration of the explicit form of the metric in global coordi- 
nates). Surfaces of constant coordinate time t are drawn, 
with the circulating arrows denoting the direction of in- 
creasing t. We will consider region I to be the "causal 
patch" of a hypothetical observer (i.e., the region lying in 
both the causal past and causal future of the observer's 
world line) in what follows. 

— 1/2 

For 3M = A + , there are also three real roots: a 
double positive root rn and a negative r nog , given by: 
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This mass is known as the Nariai mass, and in this space- 
time there is only one horizon at the positive root. The 
conformal diagram for this spacetime is shown in Fig. |3 
[l9|. There is also a time- reverse solution, starting at 

— 1/2 

past null infinity and ending at r — 0. For 3 M > A + , 
there is one real negative root, and therefore no horizons 
in the spacetime. The conformal diagram for this case is 
like figure [3 but with the horizon lines excised. 

Inside the false- vacuum region the spacetime is de Sit- 
ter space, with metric 
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there are three distinct real roots of 



in the static foliation. Fig. [3] shows the conformal dia- 
gram for the de Sitter region. Again, surfaces of constant 
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FIG. 2: Conformal diagram for the Schwarzschild de Sitter 
geometry when 3M = A+ . 




FIG. 3: Conformal diagram for the de Sitter geometry. 



coordinate time t are shown, with the arrows denoting 
the direction of increasing t. We consider region III to 
be the causal patch in which our hypothetical false vac- 
uum observer resides. 



B. Equation of Motion 



SdS exterior in Einstein's equations yields an equation of 
motion for the bubble wall of 0, Q : 



KUr) + ) - K*Ar)J) = -AirarS) 
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where KUr)±) is the extrinsic curvature tensor in SdS and 
dS respectively. In the Gaussian normal coordinates, this 
takes the form: 
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Evaluating this in metric^ the 99 and <fi(j) components 
of Eg. ITUI reduce to: 



Pds - Psds = iwcrr, 



with the definitions 
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Here, a is the metric coefficient in dS or SdS. The sign of /3 
is fixed by the trajectory because dt/ dr could potentially 
be positive or negative (motion can be with or against the 
direction of increasing coordinate time indicated in Fig^ 
and 03 . 

A set of dimensionless coordinates can be defined, in 
which Eg. 1121 can be written as the equation of motion of 
a particle of unit mass in a one dimensional potential |18| . 
Let: 



The bubble wall worldsheet has metric: 
ds 2 = -dr 2 + r(T) 2 dn 2 , 



(7) 



where r is the proper time in the frame of the wall, and 
(0, cj>) are the usual angular variables. 

The coordinates in the full AD spacetime are chosen to 
be Gaussian normal coordinates constructed in the neigh- 
borhood of the bubble wall worldsheet. Three of the co- 
ordinates are (r, 9, (f>) on the worldsheet, and the fourth, 
r], is defined as the proper distance along a geodesic nor- 
mal to the bubble worldsheet, with r\ increasing in the 
direction of SdS (true vacuum). 

The transformation from the static coordinate systems 
of SdS and dS to the Gaussian normal system can be 
constructed in closed form using the methods of 

2], 

and 

the full metric takes the form: 



ds 2 = g T r(T, rj)dT 2 + drj 2 + R(t, r)) 2 dQ 2 



(8) 



where 77 = defines the wall and therefore g TT (r, 0) = — 1 
and R(t, 0) = r(r). 

The energy momentum tensor on the wall is: 



wall 



= -a^ v 5(ri) 
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where is the metric on the worldsheet of the wall for 
\i = v = T,9,<fi and zero otherwise, and a is the energy 
density of the wall. 

Using the metric|S]and the energy-momentum tensor[5] 
together with the contributions from the dS interior and 
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where M is the mass appearing in the SdS metric coeffi- 
cient, and 
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With these definitions, Eo. 1121 becomes 
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FIG. 4: The potential for various Y. 



C. Allowed Solutions 

A bubble wall trajectory is characterized by Q =const., 
and there are three general types: 

• Bound solutions with Q < V m & x . These solutions 
start at z = 0, bounce off the potential wall and 
return to z = 0. 

• Unbound solutions with Q < V max . These solutions 
start at z = oo, bounce off the potential wall and 
return to z — oo. 

• Monotonic solutions with Q > V max . These solu- 
tions start out at z = and go to z = oo, or execute 
the time-reversed motion. 



Note that a small negative Q corresponds to a large mass, 
so that even between — I < Q < the mass will vary by 
many orders of magnitude. 

The parameter space allowed by the junction condi- 
tions is characterized by the value of the cosmological 
constant inside and outside the wall. For < A+ < A_, 
we have — 1 < Y < 1. The maximum V m3 ,x of the poten- 
tial V{z) then satisfies -2 5 / 3 - 2~ 4 / 3 < V max < 0. The 
potential curves over the entire range of Y are shown in 
FigH 

The interior and exterior cosmological constants can 
be expressed in terms of k 2 as A + = Ak 2 and A_ = Bk 2 . 
With these choices, the dynamics of the bubble wall are 
entirely determined by A, B, and Q. 

Let us now discuss some realistic values for the pa- 
rameters in this theory. The interior cosmological con- 
stant (A_ = Mf/Mt) and the bubble wall surface en- 
ergy density (k — AnMf /M 3 Z ) will be set by the scale 
of inflation (Mi). The exterior cosmological constant 
(A + = Mfc/Mfy will be set by a scale M A . These yield 

A= (47r) 2 Mf B = (4tt) 2 M 2 ' (21) 

We will consider three representative energy scales, 
covering the interesting range of energy scales for infla- 
tion. For weak scale inflation (100 GeV), k ~ 4n x 10~ 51 , 
A ~ 0, and B ~ 10 32 . For an inflation scale near the 
GUT scale (10 14 GeV > M A ), we have k ~ 4tt x 10~ 15 , 
A ~ 0, and B ~ 10 7 . Near-Planck scale inflation (10 17 
GeV) yields k ~ 4tt x 10~ 6 , A ~ 0, and S ~ 63. The 
most massive bound solution (that which just reaches 
the top of the potential) is given by converting from Q 
to M using Eq. This maximal mass is very differ- 
ent in each case, ranging from an ant-mass of M mm ~ 
10 3 M pl ~ 10~ 2 grams for Mj = 10 17 GeV to an Earth- 
mass of M max ~ 10 33 M p i ~ 10 28 grams, for M z = 100 
GeV. 



The qualitative features of a generic potential can be 
shown by considering the four illustrative (but unrealis- 
tic; see above) cases of (A = 9, B = 15) shown in Fig. 13 
(A = 1, B = 6) shown in Fig. El (A = 1, B = 2) shown 
in Fig. and (A = 2.9, 5 = 3) shown in Fig. El The 
important features are: 

• As one follows a line of constant Q in Fig. |S1 El [3 
and E| every intersection with the dashed line Q s ds 
(which is obtained by solving a s d s = for Q) repre- 
sents a horizon crossing in the SdS spacetime (this 
could either represent the black hole or cosmologi- 
cal horizons). 

• Intersections with the dashed line Qds (which is 
obtained by solving ads = for Q) as one moves 
along a line of constant Q represent the crossing of 
the interior dS horizon. 

• The vertical line on the right denotes the position at 
which Pds changes sign. /?d s > if ids is decreasing 
along the bubble wall trajectory and is negative if 
ids is increasing. For there to be a /?ds sign change, 
Y in Ea,IT9lmust be in the range — 1 < Y < jig , 
which yields the condition that B > A + 3 for a 
sign change to occur. 

• The vertical dotted line on the left denotes the 
radius at which /3 s d s changes sign. Recall that 
Psds > if i s ds is increasing along the bubble wall 
trajectory, and /3 s d s < if it is decreasing. For the 
parameters in Fig.ElandEl or (one can show) when- 
ever B < 3 (A — 1), the sign change occurs to the 
right of the maximum of V. When B > 3(A — 1) 
it occurs to the left [3^| , as shown in Fig. El (note 
that this is the interesting case where A + <C A_) 
andd 

The potential diagram contains all of the information 
needed to determine the conformal structure of the al- 
lowed one-bubble spacetimes. A complete set of the qual- 
itatively different trajectories for arbitrary interior and 
exterior positive definite cosmological constants (with 
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FIG. 5: Potential for A = 9 and _B = 15. The two dashed lines 
labeled Q s d s and Q 3 represent the exterior and interior hori- 
zon crossings respectively. The vertical dotted lines denote 
the regions in which /3 s d s and /3ds are positive and negative. 
Various trajectories are noted. 



FIG. 8: Potential for A = 2.9 and B = 3. For this choice of 
parameters, the sign change in /3 s d s occurs to the right of the 
maximum in the potential and there is no /3<js sign change. 
Various trajectories are noted. 



-0.5 



-1 





x/»,i,>0JS,d,<0 

\ 


/>„><] &<0 


, N/m-n 


\ 








— ix-m-ir/m-n 




— \ 




; iv-m-n/in 


/ 


y\ rv"-m-n"/iv-m-n 








/ 

/ 






/ 

/ 

/ 






iv-i-n/m / 




\ iv M -ni-n"/iv-i-n 


/ i 
/ / 

,./,..( 






0.5 


i 

z 


1.5 2 



FIG. 6: Potential for A = 1 and B = 6. For this choice of 
parameters, the sign change in /3 s ds occurs to the left of the 
maximum in the potential. Various trajectories are noted. 
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FIG. 7: Potential for A = 1 and 5 = 2. For this choice of 
parameters, the sign change in /3 s d s occurs to the left of the 
maximum in the potential and there is no /3d s sign change. 
Various trajectories are noted. 



A + < A_) and M > are denoted in Fig. El □ and 
[S] Figures [5] and ^| display the conformal structure of 
these solutions 39] . The conformal diagrams in each row 
are matched along the bubble wall (solid line with an ar- 
row) and the physical regions are shaded. For solutions 
with qualitatively similar SdS diagrams, the various op- 
tions for the dS interior are listed. The naming scheme 
in Fig. [S] H3 and 0D is chosen to reflect the structure of 
the conformal diagrams. The first numbers are the re- 
gions of the SdS conformal diagram that the bubble wall 
passes through. The numbers following the backslash are 
the regions of the dS conformal diagram that the bubble 
wall passes through. 

For example, consider the IV-I-II/III solution (solution 

I of Fig.[5J). These start at at r — z = 0, corresponding to 
the singularity in region IV of the SdS conformal diagram 
and the r = surface in region III of the dS conformal 
diagram. /3 s( js and j3^ s are both greater than zero over 
the entire trajectory. Therefore, the wall on the SdS 
side must follow a path of increasing t s< i s , pushing it into 
region I (crossing the black hole horizon) . The wall on the 
dS half must follow a path of decreasing coordinate time, 
and thus remain in region III. The wall then reaches a 
turning point, falls through the black hole event horizon, 
and ends up back at r = z = (the singularity in region 

II of the SdS diagram, and the r = surface in region 

III of the dS diagram). The construction of the other 
diagrams in Fig. [5] and Fig. 1101 proceeds similarly. 

Based on the inward pressure gradient, the IV-I-II/III 
solution is what one might expect the junction conditions 
to yield: a sphere of false vacuum which expands and 
then contracts. Relativistic effects, however, lead to the 
qualitatively different behavior exhibited by the other so- 
lutions in Fig. Inland Fig. ^| For instance, the evolution 
of the IV-III-II/III (solution 2) solution is qualitatively 
similar to the IV-I-II/III solution, but is so massive that 
its evolution is always hidden behind the black hole event 
horizon. The IV-TII'/III-II (solutions 3 and 4) solution 
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FIG. 9: Conformal diagrams for the allowed one-bubble spacetimes. The dS and SdS diagrams are matched across the bubble 
wall (line with arrow), and the physical regions shaded. Solutions 3 and 4 share the same SdS diagram. Regions which do not 
contain anti-trapped surfaces are shaded green (light), regions which do are shaded blue (dark). 



is a bubble which has enough kinetic energy to escape 
collapse by expanding through the cosmological horizon; 
observers inside (or who travel inside from region I) the 
false vacuum region will find themselves in an inflation- 
ary universe at late times. In the time-reverse of this 
solution, a bubble implodes from infinity into the black 
hole horizon, and the interior undergoes collapse. 

Note that there are two options for the IV-I-H'/III- 
II solution. Solution 3 corresponds to the situation 
A + 3 < B < 3(A — 1), which contains both /3 s d s and 
/3ds sign changes (see Fig. 01 . This causes the bubble 
wall to accelerate in the direction of the true vacuum 
from both the interior and exterior perspectives. Solu- 
tion 4 corresponds to the situation B < A+3, which does 
not contain a /?ds sign change (see Fig. |SJ. This causes 
the bubble wall to accelerate in the direction of the false 
vacuum from the interior perspective. The reason for the 
disparity is that solution 4 exists only when the interior 
and exterior cosmological constants are similar in magni- 
tude, and so the wall tension can have greater influence 
on the dynamics. 

The remaining solutions, shown in Fig. 1101 have in- 
teriors which approach an inflationary universe at late 
times, but lie on the opposite side of the wormhole in 



the exterior SdS spacetime. To an observer in region 
III of the SdS diagram, the IV"-III-II"/IV-I-II (solutions 
5) and rV"-III-II"/IV-III-II solutions (solutions 6 and 7) 
would appear as a "sky" of false vacuum that encroaches 
from infinity, reaches a minimum radius and then ex- 
pands back out. It has been pointed out by Bousso 
that at late times (in an asymptotically flat spacetime) 
the wall trajectory according to the observer in region 
III of the S(dS) diag ram approaches that of a true vac- 
uum bubble HI HI llj. Amusingly, the SdS observer 
will think he is in a true vacuum bubble surrounded by 
a large region of false vacuum, while the dS observer will 
think she is in a false vacuum bubble surrounded by a 
large region of true vacuum. 

This symmetry between true and false vacuum bubbles 
is made manifest in the analysis of the Coleman-De Luc- 
cia |23| instanton, which describes the production of both 
true and false vacuum bubbles 01 ■ These are zero energy 
solutions, and so we should look for an M = unbound 
solution; this corresponds (via Eg. I20fl to Q — ► — oo, and 
we see from Fig.|UEl|3 and|S]that the IV" -III-H" /IV-I-II 
solution (solution 5) or the IV" -III-II" /IV-III-II solution 
(solution 6), depending on the values of A + and A_, can 
be identified as the analytically continued false vacuum 
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FIG. 10: Conformal diagrams for the allowed one-bubble spacetimes. The dS and SdS diagrams are matched across the bubble 
wall (line with arrow), and the physical regions shaded. For solutions with qualitatively similar SdS diagrams, the various 
options for the dS interior are listed. 



instanton. The radius of the bubble wall at the turning 
point is found by considering the limit as the potential 
(Eq. I18fl goes to — oo, where on the right (unbound) side 
of the potential hump the z 2 term dominates. Solving 
for r using Eq. ^] we find the radius at turnaround to 
be 



6k 
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3k 2 



4A , A_ 



-1/2 



(22) 



which agrees with the previous literature |1(t| | (see also 
[Tl|). Since the Schwarzschild mass is zero, we are now 
matching two pure de Sitter spacetimes across the bub- 
ble wall. The conformal diagram for the exterior dS re- 
gion (right) only contains the area between the vertical 
dashed lines (which are now identified as r — surfaces) 
in solution 5 and 6 of Fig. ^| The interior half (left) 
of the diagram remains unchanged. It can be seen that 
at turnaround, the bubble will be larger than both the 
interior horizon size and the (exterior) horizon size of the 
region it has replaced. 

In the IV-III-IF/iiI-II solution (solutions 12 and 13), 
the region of false vacuum surrounding the observer in 
region III of the SdS diagram would begin very small 
and then expand out of the cosmological horizon. This 
solution will also have a time-reversed version in which 



the surrounding region of false vacuum implodes. For 
the IV"-IH"-H"/IV-III-II solution (solutions 8 and 9), 
the "sky" of false vacuum would forever reside outside of 
the horizon of a region III observer. Finally, the N/III-II 
solution (solutions 10 and 11), which we will define as 
solutions with mass greater than or equal to the Nariai 
mass, will be an exploding (or imploding in the time- 
reversed solution) bubble of false vacuum centered on 
r = 0. 

In a series of papers, Farhi et. al. p|, |2| discussed the 
application of the Penrose theorem |24j to the one-bubble 
spacetimes discussed above. If the null energy condition 
(NEC) holds (as it does for the postulated energy mo- 
mentum tensor) and there exists a non-compact Cauchy 
surface (as in the full SdS spacetime), then the existence 
of a closed anti-trapped surface in the spacetime implies 
the presence of an initial singularity. Since each point 
on the conformal diagrams in Fig. and Fig. 1101 rep- 
resents a two-sphere, such an anti-trapped surface ex- 
ists if the ingoing and outgoing future-directed null rays 
both diverge. For example, the 2-sphere represented by 
point Pi shown in Fig. solution 1, is a closed anti- 
trapped surface. This can be seen by following the null 
ray (null rays are denoted by the dashed lines in Fig. |3J) 
from r = in region III of the dS diagram to P\ and not- 
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ing that r increases monotonically as P\ is approached 
(the null rays are diverging). But following the future- 
directed null rays in the opposite direction from r = 
in region IV of the SdS diagram, across the bubble wall, 
and into the false vacuum region, shows that they also 
diverge. Thus, an initial singularity is necessary for this 
solution to exist at and near P%. This spacetime also, 
however, contains regions without anti-trapped surfaces. 
Following the future-directed null rays to point P2, for 
example, we see that the ingoing rays diverge, but the 
outgoing rays converge. For examples, see solutions 1-4 
in Fig. El where regions which contain anti-trapped sur- 
faces are shaded blue (dark) and the regions which do 
not are shaded green (light). 

If we cut the IV-I-II/III solution (solution 1) in the 
expanding phase on a spacelike hypersurface at a time 
where the radius of the bubble wall satisfies r > tbh, 
then the spacetime would not necessarily contain an ini- 
tial singularity |40| . We can remove the initial singularity 
from the IV-I-IF/III-II (solution 3 and 4) solution as well 
by cutting on the same surface. Both the IV-I-II/III (so- 
lution 1) and IV-I-IF/III-II (solution 3 and 4) solutions 
are therefore classically buildable. The IV-I-IF/III-II so- 
lution (solutions 3 and 4) is the only example where it is 
possible to form an inflationary universe from classically 
buildable initial conditions, but only exists when the inte- 
rior and exterior cosmological constant are almost equal 
(B < 3(A — 1)). This solution might be of interest in un- 
derstanding transitions between nearly degenerate vacua, 
for example in the context of eternal inflation. 

Given the existence of a classically forbidden region in 
Fig. 03 El HI and |H1 one might ask if tunneling is allowed 
from one of the recollapsing solutions (solutions 1 and 
2) to one of the expanding solutions (solutions 5-9 ) on 
the other side of the potential hump. This event, shown 
in Fig. ^2 would constitute a violation of the NEC, and 
so the Penrose theorem would no longer apply to the 
antitrapped surfaces that exist after the tunnelingevent. 
Such a process is indeed apparently allowed 0,13, and 
would describe the quantum creation of an inflationary 
universe from classically buildable initial conditions (if 
the initial condition is the IV-I-II/III solution in solution 

I)- 

This is a rather strange transition, however, as the 
unbound solution is behind the wormhole in SdS (see 
Fig. . An observer in region I would see the bubble 
expand, reach its turning point, and then disappear, only 
to be replaced by a black hole. An observer inside the 
bubble would see the wall expanding away and - just 
as it is about to turn around and start collapsing - in- 
stead disappear behind the cosmological horizon. This 
observer will be inside an inflationary universe, but for- 
ever disconnected from region I. If the black hole in the 
SdS spacetime then evaporates, the baby universe will 
become completely topologically disconnected. 

Having considered one circumstance in which an NEC 
violation precipitates the creation of an inflationary uni- 
verse, one might ask if there are others. Quantum flue- 
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FIG. 11: Tunneling spacetime 

tuations of a scalar field in de Sitter space can violate 
the NEC, and so one might imagine that any of the solu- 
tions that we have discussed which allow inflation inside 
the bubble could be spontaneously created somewhere 
along their trajectories. One example of this direct pro- 
duction of baby universes is a fluctuation into one of the 
unbound solutions (solutions 5-9). Such scenarios have 
been considered in the context of the stochastic approach 
to baby universe production by Linde |25j | and in refer- 
ence to eternal inflation by Carroll and Chen [7|. 

Another example of the direct fluctuation into an infla- 
tionary universe is the thermal decay of de Sitter vacua 
discussed by Garriga and Megevand (Ref. [23] dis- 
cusses a related mechanism). This process consists of 
the fluctuation from empty de Sitter of a bubble in un- 
stable equilibrium between expansion and collapse. This 
static solution is identified as the set of spacetimes which 
sit on top of the potentials in Fig. [SJ H3 [7| and |H] The 
two possible configurations are shown in Fig. 1121 where 
it can be seen that the bubble wall can lie on either side 
of the worm hole depending on the sign of /3 s d s at the top 
of the potential. It can be shown that the "Nariai limit" 
in Ref. (2(| corresponds to B = 3 (A — 1), where the /3 sc js 
sign change occurs at the max of the potential. 
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FIG. 12: Conformal diagrams for the thermal decay of de 
Sitter vacua. The upper diagram corresponds to B < 3(A— 1), 
the lower diagram corresponds to B > 3(yl — 1). 

With these considerations, there are at least three com- 
peting channels for the formation of baby universes: di- 
rect production, the FGG mechanism, and the zero-mass 
false vacuum instanton. It would be desirable to develop 
a scheme to directly compare the relative probabilities 
of each of these processes, and this is currently being 
explored. 

Our full catalog of solutions is also interesting in re- 
gards to a recent proposal [fH ITo| that false vacuum re- 
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gions, assumed to be larger than the interior horizon, 
must at all times be larger than the exterior, true vac- 
uum, horizon. The basis of this conjecture is the condi- 
tion that the divergence of a congruence of future directed 
null geodesies (defined as 9) must satisfy 



d6_ 

dT 



< 0, 



(23) 



where T is an affine parameter, if the NEC holds for all 
T. Null rays in the dS and SdS spacetimes satisfy this 
inequality (in dS, the inequality is exactly zero), but we 
should check that the junction conditions do not violate 
it. One requirement imposed by Eg. 1231 is that the diver- 
gence of the null rays does not increase at the position of 
the wall as they go from a true vacuum region into a false 
vacuum one. Along any given null geodesic in the bub- 
ble interior or exterior, the value of r is either increasing 
or decreasing monotonically as a function of T. We can 
therefore state the condition Eq. [231 as: one cannot have 
a null ray along which dr/dT < outside the bubble and 
dr/dT > inside the bubble. Surveying the solutions in 
Fig. Eland Fig.^J| we see that this is indeed always true. 

The authors of Ref. [iH fTEj intended to demonstrate 
that if one requires the false vacuum region to be larger 
than the interior horizon size at all times (so that in- 
flation is unstoppable), it is necessarily larger than the 
exterior horizon size. Although all of the allowed one- 
bubble spacetimes satisfy the condition Eq. [231 there are 
only two examples in which all observers agree that this 
requirement is met: the false vacuum instanton (IV- 
III-II" /IV-I-II solution (solution 5), with M = 0) and 
the IV-I-II'/III-II solution (solution 3 and 4, after turn- 
around). In every other case (including the FGG space- 
time in Fig.rrHl. the observers in region I of the SdS con- 
formal diagram will see only a black hole horizon sized 
volume replaced by the false vacuum bubble. We there- 
fore conjecture that if one requires the false vacuum re- 
gion to be larger than the interior horizon size at all 
times, then it will replace a volume larger than the exte- 
rior horizon size according to only some observers. If one 
relaxes this requirement, then there are a diverse range of 
solutions which might describe the spawning of an infla- 
tionary universe. For example, the bubbles in solutions 
3, 4 and 10-13 of Fig. l9l and 1101 grow from an arbitrarily 
small size. 

Having discussed the character of the various solutions, 
we turn now to a potentially dangerous detail, which is 
particularly important for the FGG mechanism: there 
exists a classical instability against aspherical perturba- 
tions in the spherically symmetric solutions to the junc- 
tion conditions. 



III. PERTURBATIONS 

The solutions described in Sec. Ill CI assume that the 
region of false vacuum is spherically symmetric. The sta- 
bility of these solutions against aspherical perturbations 



has important consequences, especially if one hopes to 
build plausible cosmologies. That there might be an in- 
stability in domain walls was first discussed by Adams, 
Freese and Widrow |2^. The bubble wall can trade vol- 
ume energy for surface energy and wall kinetic energy 
locally as well as globally, and so the bubble wall will 
become distorted if different sections of the wall have 
different kinetic energies. As long as the local distor- 
tions of the wall remain small compared to the size of 
the background solution's radius, this process can be for- 
mulated quantitatively as perturbation theory around a 
background spherically symmetric solution. 

Previous authors [23,|2H3(| have considered perturba- 
tions on expanding bubbles of true- vacuum [2l|, |22, |23| , 
which have zero total energy (surface, volume, and ki- 
netic energies canceling), and so can expand asymp- 
totically. As was first pointed out by Garriga and 
Vilcnkin [2£j, even though local observers on the bub- 
ble wall see perturbations grow, external observers see 
them freeze out because they do not grow faster than 
bubble radius. 

The story is different for the bound (solutions 1 and 
2) and unbound (solutions 5-9) false vacuum bubbles: 
since they reach a turning point, the perturbations have a 
chance to catch up to the bubble's expansion and become 
nonlinear. This also presumably has implications for the 
thermal decay mechanism of Garriga and Megevand [2(| , 
depending on the duration of time the bubble wall sits in 
unstable equilibrium between expansion and collapse (see 
discussion in Sec. Ill C|) . The remainder of this work will 
focus on the instability of the bound IV-I-II/III solution 
(solution 1), since physically plausible initial conditions 
may be clearly formulated. There is no obvious set of 
initial conditions for the perturbations on the unbound 
solutions, and so we simply observe that the results we 
will obtain for the bound solutions apply qualitatively 
here as well. 

To simplify the problem, we assume that the full grav- 
itational problem described in the previous sections can 
be treated as motion of the bubble wall in a fixed SdS 
background. This assumption must be validated (as we 
do below), but we are mainly interested in the low-mass 
bound solutions for which we might expect the gravita- 
tional contributions to be small. Assuming that a thin 
spherically symmetric bubble wall separates an internal 
dS from an external SdS spacetime, we can employ the 
action El El Hi: 



S = -c 



d^Xs 



(24) 



where a is the surface energy density on the bubble wall, 
Jab (a, b = 1,2, 3) is the metric on the worldsheet of the 
bubble wall, e is the difference in volume energy density 
on either side of the bubble wall: 



A+ - A_ 

and g a /3 is the metric of the background spacetime. 



(25) 
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A. Wall Equation of Motion 



The equation of motion resulting from Eq. [2] is : 



„ a b is 

a 



(26) 



where K ab is the extrinsic curvature tensor of the world- 
sheet of the bubble wall, 



K a b = -d a x^d b x v A,ra M , 
and D v is the covariant derivative and n„ 



(27) 

is the unit 
normal to the bubble wall worldsheet. 

We will use the static foliation of the SdS spacetime 
(see Eq. ^| as the coordinates x^ for the background 
spacetime. The world sheet is given coordinates (r, 9, ct>) 
as in Eq. |7| and has metric: 

lab = g^d a x f *d b x 1 ', (28) 

with the gauge freedom in choosing t fixed by 



dt =+>- 
dr~ 



V a + r' 2 



(29) 



so that 7 rr = — 1. Here and henceforth primes will de- 
note derivatives with respect to r. The other non-zero 



components of 7 a t> are 'jeg 



and 



r 2 sin 2 i 



The first task at hand is to find the worldsheet 's unit 
normal, which by spherical symmetry has only r and t 
components. Requiring orthogonality to the worldsheet 
{g^vn^daX^ = 0) and unit norm {g^n^n" = 1) yields its 
components: 



n t — —r', n r = t . 
The components of K ab are given by 



(a + r' 2 )- 1 / 2 , 



„ 1 da 
+ 2d^ 

rat' sin 2 6 = Kao sin 2 



(30a) 



(31a) 



(31b) 



Substituting Eq. into Eq. |2f)] gives the equation of 
motion for the bubble wall: 



1 da 



S' = ly/c^-±(a + t*)-~. (32) 
a r 2 dr 



Eq. ^| supplies the velocity of the bubble at some po- 
sition along its trajectory 



z =[Q-V(z )} 1 / 2 . 



(33) 



Choosing this boundary condition is effectively restrict- 
ing ourselves to the IV-I-II/III (solution 1) or IV-III- 
II/III (solution 2) solutions. Since the solutions to 
Eq. 1321 approximate the dynamics of the junction con- 
dition problem, we should parametrize by A, B, and Q. 



This can be done by using the conversions defined in 
Sec. I11BI and gives: 



z = — - 



3(5 - A) 



v / a(-Q) + z' 2 



(34) 



2,, n ., m (-Q)da 

-H-Q) + z )-— 



where a is written in terms of z as 



a = 1 - 



12 



12A 



cz(-Q) c 2 (-Q) ' 



and 



c= [\(A + B + 3) 2 -4AB\Y 



(35) 



(36) 



To justify the use of the simplified dynamics described 
above, Eq. 1341 was numerically integrated, and the po- 
sition of the turning point compared to the correspond- 
ing point on the full junction condition potential. Over 
the range of Q corresponding to the bound solutions, we 
find excellent quantitative agreement (well within 1%) 
between the turning points of the solutions to Eq. 1541 
and the junction condition potential. This was repeated 
with equally good results for the weak, GUT, and Planck 
scale potentials and also for various initial positions be- 
tween the black hole radius and the potential wall (turn- 
ing point). This shows that to zeroth order, dynamics as 
motion in a background is valid, and strongly suggests 
that it will be at higher orders as well. 



B. Perturbations 

We are now in a position to discuss the first-order per- 
turbations on the spherically-symmetric background so- 
lutions discussed in Sec. IIII Al Physical perturbations 
are normal to the worldsheet of the (background) bubble 
wall, and can be described by scalar field 4>{x) by taking 
the position of the perturbed worldsheet to be 



(37) 



where x M is the spherically symmetric solution, and 
is the unit normal to the worldsheet. It is assumed that 
cj) is much smaller than the bubble wall radius, so that a 
perturbative analysis can be made. 

The equation of motion for the perturbation field <p(x) 
in a curved spacetime background can be derived from 
the action Eq. [Rafter expanding to second order in 4>(x) 

M 



Ac 



where 



-R^hf + i? (3) - ^ 0-0, (38) 



A0 = —d a (V=ir b d b cj>) , (39) 
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and W is 



The components of the Ricci tensor are given by: 



(40) 



To solve the equation of motion, we can decompose 
(x) into spherical harmonics 



(41) 



and separate variables to get an equation for <pi m (r) . The 
geometrical factors in Eq. |2H1 become dependent on 9 or 
<f> only at second order, so we will always be able to make 
this decomposition. 

is then given by: 



2r> 1(1 + 1) 



'Jim- 



The components of h/*" are: 



h" 
h 66 



a + r 



12 



(42) 

(43a) 
(43b) 



a ^2 a r\ 1 

Rtt — -^o r a H — o r a = —a~*R rr: 



R</>c/> = Res sin 2 6 = (1 — a — rd r a). 
Contracting equations 1431 and 1441 gives: 



fi ^ = 2(l-a)_3^a_g 2 a =3 
r z r 2 



The Ricci scalar on the world sheet is 



i?(3) = 4(i + r '2 + 2 rr"), 



(44a) 



(44b) 



(45) 



(46) 



where r" is given by Eq. 1321 

After substituting Eq. E3 Eq. gSJ and Eq. gH| into 
Ea . 1381 the equation of motion for ^ m (r) is 



r e 2 



4e 



(a + r' 2 ) 1/2 + 3A + + ^ + 



2 da 6r' 2 2(1 - 4a) /(/ + 1) 



(T z c?' r ar r z 

In terms of the dimensionless variables of the junction condition problem this reads: 
2z' 



2r'c 



108A 9(A-B) 2 12(B - A) 



+ 



cz 



(a(-Q) + z' 2 ) 1/Z + 



/a , 2(-Q) 



(47) 



(48) 



(4a -1) + 



6z' 2 2(-Q)da Z(Z + 1)(-Q) 



where $ is the dimensionless perturbation field defined 
similarly to z (see Eq. I14f) . The first term acts as a 
(anti)drag on (shrinking) growing perturbations. The 
last term in this equation is always negative, acting as a 
restoring force. Perturbations will grow when the other 
terms (which are positive over most of the trajectory in 
the expanding phase) in this equation dominate. Further, 
the last term indicates that lower I modes will experi- 
ence the largest growth. The full details of the solutions, 
however, require a numerical approach, to which we now 
turn. 



been investigated only under the assumption of spherical 
symmetry, and this would be grossly violated if pertur- 
bations on the bubble wall become nonlinear. In this sec- 
tion, we investigate the circumstances under which this 
is the case. The two basic questions at issue are: first, 
when do perturbations go nonlinear for some given set of 
initial perturbations, and second, what initial perturba- 
tions can be expected. 



A. Dynamics of the Perturbation Field 



IV. APPLICATION TO THE 
FARHI-GUTH-GUVEN MECHANISM 

The possibility of creating an inflating false- vacuum re- 
gion via quantum tunneling (described in Sec. Ill C|l has 



Let us begin with the first issue. Since Eq. HHj is a 
second order ODE, it can be decomposed into the sum 
of two linearly independent solutions 

$im(T) = $ ]m (T = 0)f(l,z o ,Q,T) 

+& lm (T = 0)g(l,z ,Q,T). (49) 
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The functions f(l,Zo,Q,T) and g(l,z ,Q,T) can be 
found numerically by alternately setting $i m (T = 0) and 
$[ m (T = 0) to zero, then evolving the coupled Ea.l48land 
Eq. 1341 numerically for a time T with initial conditions 
for Q, zq, and I. If the bubble is to tunnel, it will do so 
at the time T max , when the bubble reaches its maximum 
radius and begins to re-collapse. Given / and g at time 
Tmax, the size of the perturbations at the turning point 
for any zq, Q, I, $im(T = 0), and <f>( m (T = 0) can be 
determined. An RK4 algorithm with adaptive step size 
was used to solve for / and g, with numerical errors well 
within the 1% level. 

The results of this analysis for I — 1 and for the 
low (weak) and intermediate (GUT) inflation scales dis- 
cussed below Eq. [2] are shown in Fig. ^| The solid 
lines show contours of constant (log) amplification factor 
/ (left) and g (right) versus the bubble starting radius z 
and mass parameter Q, with bubble mass increasing to- 
ward the top. The shaded regions indicate regions which 
we have disallowed as bubble starting radii because the 
bubble would not be classically buildable for r < tbh 
(marked as Q > Qbh), or the bubble is in the forbidden 
region Q < V(z) of the effective ID equation of motion 
Eq. El or the bubble would be too small to be treated 
classically. We choose the latter radius as fifty times the 
Compton wavelength z com pton of a piece of the bubble 
wall [41j. The choice of fifty Compton wavelengths is 
rather arbitrary; the effect of a larger bound would be 
to exclude more of the parameter space in Fig. ED This 
(unshaded) parameter space includes all classical initial 
conditions which could be set up by the observer in region 
I of the SdS conformal diagram. 

It can be seen in Fig. ^|that the growth of the per- 
turbations is in general larger for higher-mass bubbles 
(smaller \Q\, larger — log 10 (— Q)). The lower the inflation 
scale, the closer to zero the peak in the potential function 
becomes, and the smaller \Q\ (higher mass) bubbles are 
allowed, so at low inflation scales / and g can be very 
large. Growth for the Planck-scale inflation bubbles is 
very small, with / of order 10 and g of order 1, and is 
not plotted. 

The enhanced growth at small |<5| is due to the sup- 
pression of the term in Eq. proportional to 1(1 + 
1)(— Q), which always acts to stabilize the perturbations. 
Another consequence of this suppression is that the range 
in I over which solutions are unstable depends on Q; 
as a general rule of thumb, approximately a few times 
(— Q) -1 / 2 /-modes are unstable (note that this is unlike 
the case true vacuum bubbles, for which only the / = 0, 1 
modes are unstable). An example of the / function for 
Q = — 10~ 4 with the intermediate (GUT) inflation scale 
parameters is shown in Fig. The / functions for very 
large I modes are stable and approach sinusoids with am- 
plitudes less than one (see the inset of Fig. 114)) . meaning 
that the perturbations are never larger than their initial 



Initial Conditions and Evolution to the Turning 
Point 



Having fully characterized the growth of the perturba- 
tions, we now require an estimate for their initial values 
when the bubble is formed. There is no reason to expect 
that a region of false vacuum will fluctuate into existence 
with anything near spherical symmetry, nor is it likely 
to have thin walls (there is no instanton or other mecha- 
nism to enforce these symmetries). Since low-/ (relative 
to (— Q) -1 / 2 ) modes are unstable, an initially asphcri- 
cal bubble will only become more aspherical; this is in 
marked contrast to true vacuum bubbles, which both 
start spherical, and tend to become more spherical as 
they expand. 

Suppose, however, that we consider the best-case sce- 
nario in which a bubble is, by chance or design, spher- 
ically symmetric. It will nevertheless inevitably be 
dressed with zero-point quantum fluctuations of the per- 
turbation field. We may then check whether these fluctu- 
ations alone, considered as initial values for the perturba- 
tions of a bubble starting with a given Q and zq, suffice to 
make the bubble nonlinearly aspherical by turnaround. 

We assume that the ensemble average of the quantum 
fluctuations at the time of nucleation is zero; but the en- 
semble average of the square of the field (the space-like 
two-point function ($(0, 4>)$(d, <j>)) = ($$)) will not gen- 
erally vanish. We can write the mode functions (Eq. I41|) 
in terms of it as: 



(*L) 



/ 



(50) 



<mdn($$)Y lm (e,<f>)Y^(e,< 



By spherical symmetry, the two-point function can be 
written as a function of the angular separation "J" be- 
tween (0,0) and (0,</>), and decomposed into Legendre 
polynomials: 



= £ C, (cob ¥). 



(51) 



Using the addition theorem for spherical harmonics, we 
can write this as 

= £ 2j77i C ''^'m'(M)W^)- (52) 

l',m' 

Substituting this into Eq. 1501 and using the orthogonality 
of the spherical harmonics yields the relation: 



$f ) = 



4ttC, 
21 + 1' 



(53) 



Given some space-like two point function at the time 
the bubble is nucleated, we can obtain the C/ from 

21 + 1 f 1 

Ci = — / d(cos*)($$)P/(cos*) (54) 

47T 
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and therefore set the typical initial amplitudes of the 
mode functions as the r.m.s. value (^j^) 1 / 2 from Ea.1531 
The velocity field can be decomposed into spherical har- 
monics just as $ was, and the analysis performed above 



carries over exactly. The typical initial size of the velocity 
mode functions is then given by 



inAi 
21 + 1' 



(55) 



with 



-4, 



21 + 1 

47T 



d(cos*)($'$') p i(cos*). (56) 



The initial amplitudes in Ea.l53land Eq.___|can now be 

evolved to the turning point, and the mode functions re- 
summed. The ensemble average of the r.m.s. fluctuations 
in $ at any time at a given point will then be: 



(*(ir 



4tt 



(57) 



£ [Ct /2 f(l, zo, 1, T) + A] /2 g(l, z , q, T) 
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which can be evaluated at T = T max . 

A full model of the two-point functions ($<!>) and 
($'$') would involve quantizing the mode functions on 
the curved spacetime of the bubble wall worldsheet, 
which has a metric depending on z(T). Further, to treat 
large fluctuations, we would need to include non- linear 
terms in the equation of motion. The exact two-point 
function is therefore a rather formidable object to com- 
pute. As a simplified model, we will employ the two-point 
functions of a massless scalar field in flat spacetime, and 
replace the spatial distance r with the distance along the 
bubble wall ro^/. This massless scalar corresponds to the 
perturbations on a flat wall separating domains of equal 
energy density in Minkowski space [31| . Corrections to 
this picture in the presence of curvature should be small 
over small regions of the bubble wall. We are also ne- 
glecting the large difference in energy densities across the 
bubble wall, which will give the field a (negative) mass to 
first order. The apparent divergence of the correlator due 
to this negative mass will be rendered finite by the non- 
linear terms which must be introduced to discuss large 
fluctuations. In light of all these difficulties, and several 
more approximations we will make below, this should be 
considered as a first, rough estimate of the amplitude of 
the quantum fluctuations on the bubble at the time of 
nucleation. 

The space-like two-point function in Minkowski space 
at large separations is given by 



(0(x)^(y)) 



47rr ' 



(58) 



where r = |x— y|. As in the work of Garriga and Vilenkin 
|3l|. we introduce a smeared field operator to obtain a 
well-defined answer at close separations 



/ dVxy)- 

y-x|<s 



(59) 



where s is a smearing length, chosen to be the Compton 
wavelength s = <r -1 / 3 of a piece of wall, which is a physi- 
cally reasonable lower bound on the size of a measurable 
region. The smeared correlator will then be given by 



1 



,(x)&(y)) = -fj / d 2 z 

7T S ./|z_ x |< s 



(60) 



q-y|<s 



d 2 g(<M Z )0(q)), 



which evaluates at x = y to 



1 



-2/3 



(61) 



where Q = .916... is Catalin's constant. We have been 
unable to integrate Eq. E01 to obtain the exact form of 
the smeared two-point function for all r. However, it 
must smoothly interpolate between the value at r = in 
Eq. 1611 to the functional form at r > s given by Eq. 1581 



We therefore employ the following 'toy model' smeared 
field correlator 



(0(x)^(y)) 



a 



(2tt + AttQ - l)e 



-r 2 /2s 2 



(62) 



r/s + 1 



which has the correct asymptotics. With r = ro^, the 
dimensionless form of the toy correlator is given by: 



(-Q) 

(24tt) 2 



(2tt + 4tt£ 
1 



l) e -(* 2 /2)* 2 (63) 



m + i 



where R = tq/s. 

The velocity-velocity correlator can be calculated us- 
ing the Hamiltonian approach. For spacelike separations, 
this is given by 



6W(y)> = o- 



\p\d 2 p 
2(2tt) 2 



exp[ip(x-y)] (64) 



where we have introduced a hard momentum cut-off A 
to obtain a finite answer. This cutoff will correspond 
to the inverse smearing length, with higher momentum 
scales higher accounted for in the smeared operator. The 
integral can be evaluated in terms of generalized hyper- 
geometric functions as 



(^W(y)> = ^^ 2 (|i,| 



AV 



(65) 



At close separations, we construct a smeared opera- 
tor cj)' B . The expectation value of this operator at zero 
separation is 



log(64) + 2 7 
2tt 2 



(66) 



where 7 = .577... is the Euler-Mascheroni constant. 
Smoothly connecting the small r (Eq. I66|) and large r 
(Eq. I65fl behavior as in Eq. the $' (defined similarly 
to z') correlator on the bubble at the time of nucleation 
is: 



(-Q) 



12tt 



3 5 
lF2[ 2 ;1 '2 5 



R 2 ^ 



(67) 



12 7 -61og(64) -24 



-(i? 2 /2)* 2 



The integrals Eq. [12] and Eq. EHI for the correlators 
Eq. H3] and Eq. must be evaluated numerically. Cal- 
culation of the coefficients for every I and R is unfortu- 
nately unfeasible because of the highly oscillatory behav- 
ior of the integrands and sheer number of mode functions 
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that must be considered. However, we have been able to 
deduce sufficiently good approximate fits for C; and A; 
as a function of both I and R. In both cases, the power 
is dominated by a peak at I ~ R. 
The Ci are nicely fit by the function 



Ci 



(-Q) 83 y/21 + 1 
(24tt) 2 100 V2i? + 1 



(68) 



(V21 



l) 2 - (V2R+ l) 2 



The proposed fit for the Ai consists of two power laws 
matched at the R = I peak. For I < R, the best fit 
is A t = (-Q)/ 5 / 4 /(10i? 9/4 ) and for I > R, the fit is 
Ai = 3(-Q)R ls /(200l %6 ). Because these power law in- 
dices are slightly uncertain, we only count the modes with 
I < R in the Ai. This is conservative, and also justified 
because these modes will not contribute significantly to 
the sum in Eq. 1571 

With these initial conditions, we can now evolve each 
mode function using Eq. 1491 and then re-sum in Eq. 1571 
to find the average size of the fluctuations at the turning 
point. We have calculated / and g up to the / correspond- 
ing to the last unstable mode of the smallest \Q\, for all 
three inflation scales (Mj = 100 GeV, Mj = 10 14 GeV, 
and Mi = 10 17 GeV). The results for weak- and GUT- 
scale inflation arc shown in Fig. 1131 where the dotted 
line indicates the boundary of the region over which the 
perturbations become non-linear (non-linear to the left 
of the line). It can be seen that in this model, even just 
quantum perturbations of the bubble wall grow nonlinear 
in bubbles that start at radii less than about one-tenth of 
the turnaround radius; this grossly violates the assump- 
tion of spherical symmetry used in tunneling calculations. 
On the other hand, none of the parameter space in the 
high inflation scale case went non-linear, and at all scales 
there is always a region of initial bubble radii near the 
turnaround radius, for which nonlinearity never occurs. 



C. Thick Walls and Radiation 

Just as we have no reason to expect a fluctuated region 
to be spherically symmetric, we have no reason to assume 
that it will have thin walls. An analysis of thick-walled 
true vacuum bubbles was undertaken in Ref. |3ll l32l |. 
where it was found that the instabilities found in thin- 
walled case are still present in the form of deformations 
normal to the bubble profile. In the case of small false 
vacuum bubbles, there is no obvious consideration (such 
as a corresponding instanton) to supply the profile of the 
bubble wall, and so we can merely conjecture by prece- 
dent that the instability would be retained in the thick- 
walled case as well. 

Another consideration, applying to bubbles smaller 
than the false vacuum horizon size, is whether inflation is 
spoiled by non-vacuum contributions to the energy den- 
sity. The perturbations on the bubble wall translate into 



gravitational waves [33t |3J], and since the bound bub- 
ble solutions remain relatively close to their gravitational 
radius and become distorted over many different length 
scales on a relatively short time scale (see the quasi- 
exponential growth in Fig. I14|) . they will be emitters of 
copious gravitational radiation. Another problem arises 
if the kinetic and gradient energy of the field becomes ap- 
preciable in the bubble interior, either from intrusion of 
the wall (for example, imagine a bubble being pinched in 
half by some non-linear perturbation), or from particle 
production or other scalar modes propagating in from 
the wall. If the emission of energy into the interior of 
the bubble from any combination of these modes makes 
a significant contribution to the equation of state, then 
inflation will not occur. 



V. SUMMARY & DISCUSSION 

In this paper we have examined the feasibility of pro- 
ducing inflationary universes from a spacetime with a 
small cosmological constant. Reviewing the solutions al- 
lowed by the junction conditions, there are several dis- 
tinct possibilities if one allows for violations of the null 
energy condition. The first is the Farhi, Guth & Gu- 
ven (FGG) mechanism, in which (referring to Fig. [5] and 
Fig. EH) the IV-I-II/III (solution 1) or IV-III-II/III (solu- 
tion 2) solution is fluctuated in the expanding phase and 
then tunnels to an unbound solution, as shown in Fig. II II 
To the outside observer in region I of the SdS conformal 
diagram it appears that the bubble has disappeared be- 
hind the horizon, but on the other side of the wormhole, 
both an inflationary universe and a non-inflating universe 
(an asymptotically true-vacuum de Sitter region) have 
come into existence. In one case, shown in solution 5 of 
Fig. 1101 we can imagine an observer in region III, then 
take the Schwarzschild mass m — > limit (thus removing 
regions I, II and IV entirely) and interpret the space- 
time as that of the analytically continued Coleman-De 
Luccia instanton. Since the very same spacetime takes 
part in both the FGG mechanism and the Coleman-De 
Luccia false vacuum instanton, there may be some way 
to smoothly interpolate between these two processes; we 
intend to explore this idea further in future work. 

In addition to the FGG and regular instanton mech- 
anisms, there are two more possibilities. First, an infla- 
tionary universe might be directly produced by some null 
energy condition (NEC) violating fluctuation into one of 
the unbound or monotonic solutions shown in Solutions 
5-13 of Fig. ^3 To the outside observer in region I of 
the SdS conformal diagram, these solutions would be in- 
distinguishable from the fluctuation of a black hole, but 
they would secretly entail the creation of everything on 
the other side of the wormhole, as in the FGG tunneling. 
A second method of direct production is the fluctuation 
(which does not require an NEC violation) of the IV- 
I-IF/III-II solution (solution 3 and 4). These solutions 
only exist in the case where the interior and exterior cos- 
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mological constants are comparable, and so would not 
correspond to the nucleation of inflation from a universe 
like ours, but they might be of interest in understanding 
transitions between nearly degenerate vacua. 

Examining these classical solutions to first order, we 
have shown that an instability to aspherical perturba- 
tions exists in those solutions which possess a turning 
point. This includes the bound (solutions 1 and 2) and 
unbound solutions (solutions 5-9). In the latter case 
there is no clear way to set an initial radius or initial 
perturbation amplitude, so we can say only that collaps- 
ing bubbles are violently unstable. The bound solutions 
are amenable to quantitative investigation, and we have 
focused on the growth of perturbations in the expanding 
phase that precedes tunneling in the FGG mechanism. 
For bound expanding bubbles formed by the fluctuations 
of a scalar field in de Sitter space, there is no instanton to 
enforce spherical symmetry, so we would expect initial as- 
pherical perturbations to be relatively large. Since there 
is no detailed model for the fluctuating scalar field to see 
how large, we have instead calculated an estimate of the 
minimal deviations from spherical symmetry in light of 
quantum fluctuations, and present this as an extremely 
rare, best case scenario for spherical symmetry. These 
minimal fluctuations were then evolved to the turning 
point of the bound solutions, which is the point in the 
FGG mechanism where there is a chance for tunneling 
to an inflationary universe to occur. Of the three repre- 
sentative energy scales for inflation (false vacuum energy 
densities) we have studied, the evolved minimal pertur- 
bations on a Plank scale bubble remain small over most 
of the allowed parameter space, while the perturbations 
on GUT and weak scale bubbles can grow nonlinear if 
they start at a sufficiently small (< 10%) fraction of the 
turnaround radius. Thus even in the best-case scenario 
some bubbles become nonlinear, but on the other hand 



there will always in principle be some that do not. 

The instability introduces complications into the use 
of the FGG mechanism as a means of baby universe pro- 
duction. The existing calculations of the tunneling rate 
rely heavily on the assumption of spherical symmetry. 
It is unclear how to perform a similar calculation for a 
(possible non-linearly) perturbed bubble, as the number 
of degrees of freedom has drastically increased and the 
assumption of a minisuperspace of spherically symmetric 
metrics is no longer good. Further, the bubble interior 
will become filled with scalar gradient and kinetic energy 
and gravity waves, possibly upsetting the interior suffi- 
ciently to prevent inflation. One might argue that in an 
eternal universe there is plenty of time to wait around for 
a fluctuation which is sufficiently spherical. However, to 
fully understand the importance of the FGG mechanism, 
one must both have a model of the scalar field fluctu- 
ations, which would predict the distribution of bubble 
shapes and masses, and also a model for tunneling in 
the presence of asphericities. In addition, one must un- 
derstand the competition among the various processes 
which result in the formation of inflating false vacuum 
regions, and thus show that the FGG mechanism is not a 
relative rarity. These problems represent significant cal- 
culations in a theory of quantum gravity which we do not 
yet possess, but progress in these areas would undoubt- 
edly improve our understanding of the initial conditions 
for inflation. 
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